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Z.  ZMTBODUCTZOM 

The  non-intrusive  measurement  of  secondary  flows  in 
annular  turbines  can  provide  insight  into  the  design  and 
performance  of  turbomachinery.  Turbine  stage  efficiency  can 
be  improved  with  more  precise  understanding  of  the  flow 
mechanics  and  the  ability  to  numerically  predict  the  response 
to  blade  redesign  without  manufacturing  and  testing. 

Goldman  and  Seasholtz  [Refs.  1-3]  have  conducted 
continuing  research  in  the  field  of  laser-Doppler  velocimetry 
(LDV)  measurement  and  numerical  prediction  in  an  annular 
turbine  cascade.  Their  efforts  have  concentrated  on  obtaining 
comparisons  ahead,  in,  and  downstream  of  the  blade  passage. 
In  order  to  obtain  their  results,  an  extensive  optical  system 
with  correction  for  curved  surfaces  has  been  developed.  To 
obtain  two-dimensional  velocity,  they  have  measured  one 
component  ard  then  rotated  the  optics  to  resolve  the  second 
component . 

This  report  presents  the  results  of  laser-Doppler 
velocimetry  techniques  and  three-hole  pressure  probe 
measurements  in  the  wake  of  a  high  turning  turbine  nozzle. 
Only  radial  fiber-optic  probe  traverses  were  performed. 
Circumferential  surveys  were  obtained  by  manually  rotating  the 
annular  cascade  blade  row  upstream  of  the  measuring  location. 
Pressure  measurements  were  compared  to  numerical  results 


1 


obtained  using  a  three-dimensional  viscous  computer  program. 
Laser-velociaetry  techniques  using  two-component  fiber-optic 
equipment  were  developed  for  the  small  annulus  cascade 
including  two-dimensional  measurements  of  the  endwall  flow 
through  a  small  access  hole  in  the  outer  casing. 
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B.  TB8T  FACZLZTY 

Airflow  for  the  annular  turbine  cascade  experiment  was 
provided  by  a  VA-312  Allis-Chalmers  12  stage  axial-flow 
compressor  housed  in  the  Naval  Postgraduate  School 
Turbopropulsion  Laboratory.  The  compressor  was  operated  at 
12,000  rpm  at  various  discharge  pressures.  This  air  supply 
was  also  used  to  generate  up  to  Mach  0.795  through  a  107.95  mm 
(4.25  in)  diameter  nozzle  for  calibration  of  a  three-hole 
pressure  probe. 

Upstream  air  pressure  was  metered  through  the  use  of  two 
dump  valves  and  one  cutoff  valve  prior  to  a  plenum  chamber. 
Air  from  the  plenum  was  directed  to  the  test  section  through 
flow  straighteners  in  a  254  millimeter  (10  in)  steel  pipe  as 
shown  in  Figure  1.  Flow  stagnation  pressure  was  measured 
upstream  at  two  locations.  One  combination  probe  was 
connected  to  a  mercury  manometer  board  located  near  the  dump 
valve  controls.  This  pressure-setting  probe  also  provided  a 
digital  readout  of  flow  stagnation  temperature.  The  second 
probe,  and  the  static  ports  in  the  air  supply  pipe,  were 
connected  to  a  Scanivalve  for  automated  pressure  acquisition. 

LDV  seed  material,  atomized-glycerin  particles,  was 
introduced  through  an  eight  millimeter  (5/16  in)  diameter 


I 

I 

I  copper  pipe  as  shown  in  Figure  1.  Seed  atomization  was 

j  accomplished  using  laboratory  compressed  air  with  a  commercial 

TSI,  Inc.,  Six-Jet  Atomizer.  Laser  traverse  table  buffet  from 
the  exiting  flow  was  eliminated  with  a  prefabricated  aluminum 
cone  (Figure  1) .  Flow  characteristics  were  examined  with  the 
cone  on  and  off,  using  pressure  measurements  and  dye  injection 
from  the  two  chord  inner  hub  static  ports. 


Figure  1.  Supply  Pipe  Schematic 


B.  AMMULAR  TURBINE  CASCADE 
1.  GENERAL  DESCRIPTION 


The  Annular  Turbine  Cascade  (ATC)  was  developed  to 
provide  a  small-scale  device  for  flow  experiment  and 
computation  on  a  transonic  turbine  stator.  Thomas  [Ref.  4] 
developed  the  basic  apparatus.  Figure  2  shows  a  cutaway  view 
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Figure  4.  Wake  Positioning  System 
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The  original  ATC  pressure  neasurenent  access  was 
located  two  axial  chords  downstream  of  the  blade  trailing 
edges.  Figure  5  details  the  old  and  new  probe  holes.  The 
outer  casing  was  modified  to  accommodate  the  half -chord  access 
hole  and  is  shown  in  conjunction  with  the  adjustable  holder. 
A  new  set  of  static  pressure  ports  was  added  to  the  inner  htib 
as  shown.  The  ports  were  set  90  degrees  apart  and  in  axial 
alignment  with  the  original  holes.  Tiibing  was  connected  to 
provide  an  average  pressure  tap  for  both  hub  axial  locations. 
Appendix  A  contains  a  complete  set  of  the  drawings  used  for 
machining. 


Figure  5.  Pressure  Measurement  Changes 


3.  LDV  1I00XFZCATZ01I8 


LDV  access  was  changed  to  Bininize  light  reflection 
and  iaprove  off-axis  options.  The  desire  to  ultinately  obtain 
a  third  component  necessitated  an  offset  angle  for  a  second 
fiber>optic  probe.  The  design  allowed  for  simple  removal  and 
replacement  of  the  window  during  operation  for  cleaning,  and 
radial  access  as  shown  in  Figure  6.  The  position  of  the  wake 
relative  to  the  probe  volume  was  controlled  externally  by 
rotating  the  blade  row,  inlet  bellmouth  and  center  body  ahead 
of  the  measurement  volume.  An  alignment  blank  was  used  to 
orient  the  laser  traverse  axis  to  the  ATC.  The  1.0922  mm 
(0.043  in)  diameter  hole  in  the  blank  provided  an  endwall  flow 
measurement  capability  as  shown  in  Figure  7. 
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Figur«  7.  Endwall  Measurement  Technique 

C.  PSBSSURB  PROBE 

Total  pressure,  Mach  number  and  flow  angle  measurements 
were  made  with  a  three  hole  yaw  probe  [Ref  4].  All  pressure 
taps  were  connected  to  a  Scanivalve  with  a  ±344  kPa  (±50  psl) 
differential  transducer.  Scanivalve  control  and  data 
recording  were  accomplished  using  an  HP-9000  computer  system. 
Appendix  B  contains  calibration  curves  and  a  listing  of  the 
Hewlett-Packard  code  used  to  acquire  and  output  data. 
Pressure  data  from  the  ATC  were  acquired  with  the  same 
hardware  and  software  used  in  the  probe  calibration. 

Previously,  yaw  alignment  was  accomplished  using  a  "U"- 
tube  water  manometer.  Difficulty  encountered  with  hysteresis 
and  the  response  of  the  water  to  yaw  misalignment  was 
corrected  wluh  the  use  of  a  differential  pressure  transducer. 


D.  XASm-OOPPLIB  VBLOCZlOm 


The  LDV  system  is  shown  schematically  in  Figure  8.  The 
Lexel  Model  95  four-Watt  Argon-Ion  LASER  was  connected  to  a 
TSI,  Inc.,  Model  9201  ColorBurst  Multicolor  Beam  Separator. 
The  beam  separator  split  the  incoming  light  into  two  beams, 
with  one  beam  passing  through  a  Bragg  cell  for  40  Mhz 
frequency  shifting.  The  two  beams  were  further  split  into 
three  polarized  pairs:  green  (514.5  nm) ,  blue  (488  nm) ,  and 
violet  (476.5).  The  TSI,  Inc.,  Model  9230  ColorLink 
Multicolor  Receiver  provided  the  frequency-shift  signal  to  the 
ColorBurst  for  beam  shifting.  Each  beam  was  directed  through 
Individual  couplers  for  transmission  to  the  laser  probe  by 
fiber-optic  cable.  The  fiber-optic  probes  contained  receiving 
optics  to  direct  the  return  signals  to  the  ColorLink. 
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tolorLink  provided  {diotoaultipller  and  frequency 
ehifting  functions.  Control  of  shifting  and  photomultiplier 
voltage  was  accomplished  remotely  by  computer.  After 
conditioned  by  the  ColorLink,  the  signals  were  sent  to 
individual  TSI,  Inc.,  Model  1990C  signal  processors  where 
valid  Doppler  signals  were  counted  with  reference  to  desired 
characteristics  set  by  the  user.  An  MI-990  interface 
monitored  coincidence  criteria  (2-D  measurements)  and 
transferred  valid  counter  information  to  computer  memory. 

The  fiber-optic  probes  were  mounted  to  a  Model  41583 
LINTECH  Traverse  Table.  An  Applied  Motion  Products  System 
1618  Traverse  Controller  was  used  to  control  traverse  movement 
manually  and  by  computer.  [Ref.  5]  All  remote  control  of  LDV 
system  equipment  and  data  processing  was  completed  through  the 
menu  driven  environment  of  TSI's  FIND  version  4.04. 

E.  EXPERIMENTAL  PROCEDURE 
1.  PROBE  CALIBRATION 

Probe  calibration  and  investigation  of  data 
acquisition  options  were  completed  together.  The  calibration 
was  accomplished  using  a  free-jet  with  a  107.9  mm  (4.25  in) 
nozzle  connected  to  the  air  supply  (Figure  1)  used  for  the 
annular  cascade.  Figure  9  shows  the  testing  apparatus  and 
pressures  used  to  calibrate  the  cobra  probe. 

The  probe  was  aligned  in  yaw  with  a  ±68.9  kPa  (±10 
psi)  Stathan.,  SCO  277-10,  differential  pressure  transducer. 
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Yaw-allgiment  response  was  an  improvement  over  the  previous 
manometer  arrangement.  The  transducer  bridge  output  was 
periodically  zeroed  by  removing  the  probe  from  the  flow.  The 
yaw  pressure  measurement  was  tapped  from  one  side  of  the 
transducer  as  shown  in  Figure  9.  The  Scanivalve  system 
recorded  the  five  pressures  shown  in  Figure  9.  Pre-operation 
tubing-leak  checks  were  completed  by  trapping  pressure  to  each 
Scanivalve  port  in  turn. 


The  probe  was  calibrated  at  Mach  numbers  from  0.20  to 
0.79  in  steps  of  0.1.  Pitch  sensitivity  was  investigated  at 
twelve  positions  from  minus  ten  degrees  to  plus  eight  degrees 
(Figure  9) .  Appendix  B  contains  a  description  of  the 


calibration  variables  and  shows  the  sensitivity  of  the  probe 


to  plt^  anglo.  Sixth-order  polynomials  were  fit  to  the 
calibration  data.  The  zero-degree  probe  calibration  was 
successfully  validated  (Appendix  B)  under  different 
atmospheric  conditions. 

2.  PRBS8URB  PROBB  KEABUSBMEIITB 

Wakes  were  surveyed  by  repositioning  the  blade  passage 
upstream  relative  to  the  probe.  The  inner  hub  rotated  with 
the  blade  ring  as  the  wake  was  positioned  for  LDV  or  pressure 
measurement.  The  bolt  shown  in  Figure  4  was  used  to  position 
and  hold  the  inner  ring  and  stator.  A  total  of  18  degrees 
circumferential  variation,  in  steps  of  one  degree,  was 
achieved. 

Cobra-probe  measurements  were  completed  as  described 
in  the  calibration  section.  Upstream  static  and  stagnation, 
and  inner-hub  static  pressures  were  used  in  conjunction  with 
the  center-port  and  yaw-port  pressures  of  the  probe.  The 
probe  provided  total  pressure  and  flow  angle.  Mach  number  was 
obtained  from  P  (Appendix  B) .  The  desired  hub-static  to 
upstream-stagnation  pressure  ratio  was  set  using  the 
Scanivalve  system  and  manual  control  of  the  air  supply  dump 
valves.  The  probe  was  positioned  radially  and  aligned  with 
the  flow  before  pressure  data  were  acquired.  Upstream 
stagnation  temperature  and  wake  flow  angle  were  recorded 
manually  for  each  measurement  location.  Flow  angle  was  read 
from  a  vernier  scale  to  an  accuracy  of  plus  or  minus  0.2 
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d*gr«e8  at  span  locations  greater  than  25  percent.  Span 
position  was  est2d>llshe<}  with  reference  to  100  percent  at  the 
outer  wall. 

3.  LD7  hLZOmiSlIT 

LDV  alignment  for  endwall  flow  measurement  was 
attempted  as  shown  in  Figure  10.  The  laser  mounting  bracket 
micrometer  was  set  to  ten  millimeters  for  visual  alignment 
around  the  access  hole  as  shown.  Obtaining  a  repeatable  depth 
reference  was  not  possible.  This  alignment  technique  did 
provide  micrometer  scale  coverage  for  data  recording  and 
ensured  the  beams  passed  cleanly  through  the  hole.  The  length 
of  the  probe  volume  precluded  the  use  of  the  outer  surface  of 
the  plate  as  a  reference  point. 


Figure  10.  Endwall  LDV  Alignment  Technique 
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Alignment  for  automatic  traverse  of  the  wake  passage 
through  the  window  was  accomplished  using  the  alignment  blank 
(Figure  A8) .  The  multi-component  beams  were  positioned  in 
each  of  the  six  alignment  holes  using  manual  control  of  the 
traverse  table.  The  coordinates  were  read  from  the  digital 
position  indicator  on  the  traverse  controller  for  each  hole. 
Axis  conversion  from  the  traverse  to  the  ATC  was  completed 
using  the  position  data  and  the  dimensions  of  the  alignment 
blank  (Appendix  A) .  Traverse  table  movement  could  then  be 
programmed  using  the  automatic  traverse  control  menu  in  TSI's 
FIND  software. 

4 .  ENONALL  MEASUREMENTS 

Endwall  flow  data  were  acquired  without  the  use  of 
automatic  traverse  control.  The  laser  beams  were  aligned  with 
the  half-axial-chord  access  hole  as  described  above,  and  then 
traversed  manually  using  the  probe-bracket  micrometer. 
Adjustments  to  seeding  position,  counter  gain  and 
photomultiplier  voltages  were  continually  updated  to  maximize 
data  rate  and  minimize  the  effect  of  noise.  Seed-wand 
position  was  particularly  important.  Micrometer  position 
information  was  married  to  the  appropriate  data  files  in  post 
processing  of  FIND'S  output  ASCII  .PRN  files  using  Borland 
Corporation's  spreadsheet  program  Quattro  Pro  version  5.0. 
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ZXI.  CONPUTATZOIIAL  FLUID  DYMAMICS 

A.  OAZO  OBMBRATZOM 

Grid  generation  was  completed  using  Turbomachinery  C  GRID, 
TCGRID, [Ref.  6].  The  code  was  developed  to  build  3-D  "C  or 
"H"  type  grids  for  turbomachinery  in  the  standard  PL0T3D 
format  (NASA  Ames  Research  Center  graphics  software) . 
Appendix  C  contains  the  namelist  input  file  used  to  generate 
the  stator  g.iometry  and  a  discussion  of  the  output.  The  inlet 
bellmouth  and  center  body  as  shown  in  Figure  4  were  not 
modeled  for  this  study. 

Careful  manipulation  of  the  namelist  variables  was 
required  to  obtain  a  smooth  solution.  The  grid  geometry  was 
developed  to  Improve  resolution  of  the  endwall  boundary  layer. 
The  ability  to  generate  a  fine  grid  was  hampered  severely  by 
the  length  of  wake  desired  in  the  model.  Reduction  of  the 
wake-exit  plane  from  one  axial  chord  downstream  of  the  blade 
trailing  edges  [Ref.  4]  to  one-half  axial  chord  resulted  in  a 
dramatic  increase  in  grid  resolution.  The  final  grid 
dimension  was  150x31x65,  up  from  121x31x21.  The  blade  surface 
received  an  additional  29  points  and  the  radial  resolution  was 
improved  from  21  points  to  65.  The  smaller  computational  wake 
length  also  reduced  the  interval  between  grid  lines  from  the 
trailing  edge  of  the  blade. 
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B.  OOKPOTBTZOIIAL  BCBENE 

Flow  solutions  were  completed  using  Rotor  Viscous  Code  3-D 
(RVC3D) .  RVC3D  [Ref.  7]  was  designed  to  solve  three- 
dimensional  viscous  flows  in  turbomachinery.  The  code  solves 
thin-layer  Navier-Stokes  equations  in  Cartesian  coordinates 
using  an  explicit  finite  difference  technique.  Turbulence 
effects  were  modeled  using  the  Cebeci-Smith  option.  The 
restart  option  of  the  code  was  utilized  throughout  the 
investigation  to  monitor  convergence  of  the  solution. 
Appendix  D  contains  the  namelist  input  file  used  to  obtain  the 
final  solution  and  a  description  of  steps  required  to  run  the 
code  on  the  Naval  Postgraduate  School  (NPS)  Cray  computer. 
Output  solution  files  (fort. 3)  were  examined  using  the 
graphics  software  PL0T3D.  Comparison  with  experiment  required 
data  reduction  of  the  output  file.  Appendix  E  contains  a 
synopsis  of  the  output  file  structure  with  data  reduction 
steps  required  for  comparisons  and  a  listing  of  the  FORTRAN 
program  (plane. f)  used  to  accomplish  the  manipulation. 

The  final  solution  selected  for  comparison  was  chosen  for 
the  minimum  residual  and  the  smoother  transition  of  a  solution 
across  the  wake  cut.  An  area  of  unsteady  flow  became 
prominent  near  the  h\ib  and  on  further  analysis  revealed  an 
area  of  predicted  reverse  flow.  Figure  11  contains  a 
comparison  between  total  pressure  ratio  results  at  the  50 
percent  span  location  for  1200  and  1490  iteration  solutions. 
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One  blade  wake  is  shown  with  the  wake  cut  in  the  grid  clearly 
evident  in  the  plot.  The  closely>spaced  points  in  the  wake 
cut  of  the  TC6RID  geometry  has  introduced  a  mathematical 
boundary.  The  1200  iteration  solution  exhibited  the  beginning 
of  the  reverse  flow  area  near  the  hub,  but  contained  serious 
discontinuities  because  of  the  solution  scheme.  The  1490 
iteration  total-pressure-'ratio  plot  demonstrates  the  smoothing 
that  took  place  as  the  solution  progressed. 


CFD  Solution  Comparison 

Total  Pressure  Ratio  (50  %  Span) 


Figure  11.  CFD  Solution  Comparison,  50%  Span 
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ZV.  BB81ILT8  AMD  DZ8CD88ZOM 


A.  FR888U8B  P80BB 

Annular  Turbine  Cascade  pressure  measurements  were 
obtained  at  an  average  Reynolds  number  of  650000  based  on 
axial  chord.  The  pressure  ratio,  inlet  total  to  exit  static 
at  the  hub,  was  set  at  0.68  to  enable  a  proper  comparison  to 
be  made  with  the  CFD  solution.  Pressure  and  flow  angle 
measurements  were  repeatable,  with  Scanivalve  pressure- 
acquisition  standard  deviations  of  less  than  0.20  percent. 
However,  data  obtained  inside  25  percent  span  were  at  low 
dynamic  pressure  and  too  difficult  to  resolve  with  the 
pressure  probe,  and  are  therefore  not  reported.  Flow 
vmsteadiness  increased  as  the  probe  was  immersed  further  in 
from  the  endwall.  The  direction  of  approach  to  zero  yaw  in 
the  less-than-2 5-percent  span  region  affected  the  flow  angle 
measurement  by  as  much  as  0.5  degrees. 

Pressure  and  flow-angle  data  were  taken  at  18  wake 
positions  for  25,  50,  75  and  90  percent  span  locations  on  the 
half-axial-chord  wake  measurement  plane.  Pressure  data 
repeatability  was  demonstrated  along  the  50  percent  span 
station  for  all  wake  positions,  and  radially  at  the  two-chord 
zero-degree  wake  position.  Figure  12  shows  the  total - 
pressure-ratio  data  obtained  for  the  50  percent  span  location. 


! 
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n)«  data  ahow  that  thara  vaa  a  nagllgibla  influanca  of  tha 
conical  diffuaar.  Figure  13  ahowa  tha  repeatability  obtained 
in  radial  aurvays  for  a  conatant  wake  poaition  at  the  two- 
chord  Baaaurement  location.  The  repeatability  shown  in 
Figures  12  and  13  clearly  validated  the  ability  to  reset  the 
deaired  operating  preseure  ratio  for  the  experiment.  Flow 
angle  meaeurements  were  repeated  at  the  90  percent  span 
location,  and  radially  at  the  two-chord  wake  measurement 
location  as  shown  in  Figures  14  and  15.  In  all  cases,  the 
repeatability  of  the  probe  pressure  measurements  was 
excellent,  with  the  largest  discrepancy  less  than  0.5  percent 
in  total  pressure. 
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Total  Pressure  Ratio  Repeatability 

50%Span 


Flgur*  \2,  Pt/PO  Repeatability,  50%  Span 


Total  Pressure  Ratio  Repeatability 

2  Chord  (Diffuser  Ofl) 


Figure  13.  Pt/PO  Repeatability,  Radial 
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B.  UUnm  MBMUBBIIBBTS 

Two-HliBenslonal  LDV  neasureaents  were  obtained  through  a 
1.0922  an  (0.043  in)  hole  to  a  depth  of  approxlaately  three 
ailliaeters.  Initial  aeasureaents  were  taken  with  one 
coaponent  (green)  of  the  two-coaponent  probe.  Figures  16  and 
17  show  the  velocity  profile  and  turbulence  intensity  results 
froa  the  first  atteapt  at  endwall  aeasureaents.  The  'depth' 
shown  in  the  graphs  is  an  estiaated  distance  from  the  endwall 
based  on  the  rapid  loss  of  data  rate  observed  at  the  depth 
shown  as  0.1  mm  from  the  endwall.  The  data  map  a  boundary 
layer  profile  to  approximately  0.6  millimeters.  The  high 
pressure  ratio  (0.98)  provided  near  perfect  seed  conditions, 
and  the  casing  vibrations  were  low.  The  data  shown  for  V~z  do 
not  necessarily  coincide  with  the  exact  locations  of  the  data 
for  V-0  because  the  fiber-optic  probe  was  rotated  90  degrees 
and  alignment  could  not  be  maintained. 

The  turbulence  intensity  results  are  with  respect  to  the 
maximum  velocity  at  the  furthest  measurement  location  into  the 
annulus.  The  significantly  higher  turbulence  intensity  in  the 
V-6  direction  is  probably  due  to  the  steep  gradient  of  the  V-0 
mean-flow  profile  in  the  radial  direction. 
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Vigur*  1€.  Single  Conponent  LDV  Endvall  Velocity 


I  -o-  MlNoJTwfauttne*  VWTwtiilwK*  | 


Figure  17.  single  Component  LDV  Turbulence 
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Two  coBponent  LDV  neasureaents  were  obtained  at 
approxlnately  15  percent  higher  flow  velocity.  Figure  18 
depicts  a  test  section  'sweep'  of  near-endwall  velocity 
measurements  at  a  pressure  ratio  of  0.974.  Depth  reference 
(in  mms  with  respect  to  an  unkno%m  datum  on  a  vernier  scale) 
is  shown.  The  exact  distance  from  the  wall  was  not 
obtainable.  Flow  periodicity  is  suggested  and  decreasing 
velocities  toward  the  outer  casing  show  the  endwall  boundary 
layer.  The  )>oundary  layer  is  distorted  due  to  corner  vortices 
and  secondary  flows  that  originate  in  the  blade  row.  Flow 
angle  and  turbulence  intensity  plots  are  included  in  Figures 
19  through  21.  Once  again  the  turbulence  intensity  in  the  V-6 
direction  is  higher  than  that  in  the  V-z  direction.  Figures 
18  through  23.  represent  automatic  editing  of  the  histograms  by 
the  FIND  software  (user  selection) .  Data  outside  two  standard 
deviations  were  discarded. 

Increased  flow  velocity  introduced  higher  vibration  levels 
and  temperature  increases  as  the  test  progressed.  Hole 
alignment  was  periodically  verified  and,  occasionally,  small 
adjustments  (less  than  0.01  mm)  were  required. 
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Tlgur*  18.  Two  Component  LDV  Endwall  Velocity 


LDV  Endwali  Flow  Angle 

Depth  Reference  (mm) 
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Figure  19.  Two  Component  LDV  Flow  Angle 
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Vttieta  Turbulence  Intensity 
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Figure  20.  Two  Component  LDV  Turbulence,  V-theta 


Vz  Turbulence  Intensity 
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Figure  21.  Two  Component  LDV  Turbulence,  V-z 


c.  MinfBszau.  sbsults 

The  CFD  solution  was  run  with  the  same  pressure  ratio 
(0.68)  used  in  the  experiment.  Graphical  output  of  the 
numerical  results  was  useful  in  obtaining  a  qualitative 
understanding  of  the  flow  characteristics.  Figure  22  depicts 
the  predicted  pressure  ratio  for  the  blade  surface  at  the  hub, 
mid-span  and  tip.  The  suction  surface  curves  suggest  that  the 
position  of  the  blade  passage  throat  was  at  0.8  of  axial 
chord.  Figure  23  shows,  with  colored  contours,  pressure 
coefficient  information  for  the  blade  surface  pressure 
distribution.  The  Mach  number  contours  shown  in  Figure  24 
provided  an  initial  indication  of  where  flow  characteristics 
were  well  established.  Boundary  layers  appear  as  expected, 
thin  near  the  endwall  and  thick  near  the  hub.  Wake  curvature 
in  the  spanwise  direction  is  also  evident  in  the  Mach  contour 
plot.  Difficulty  in  obtaining  repeatable  flow  angles  in  the 
experiment  below  25  percent  span  may  have  been  due  to  the 
unsteady  f].ow  separation  which  is  suggested  by  the 
computational  results  for  the  region  near  the  hub  wall. 
However,  the  reverse  flow  predicted  by  the  code  in  this  region 
could  not  be  verified  in  the  experiment. 
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Figure  23.  PL0T3D  Blade  Pressure  Coefficient  Plot 
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.gure  24.  PL0T3D  Exit  Plane  Mach  Contour  Plot 
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D.  COMPARISON 


Three  comparisons  were  made  between  numerical  and 
experimental  results.  To  show  the  flow  periodicity  ,  the  CFD 
solutions  were  repeated  over  two  blade  passages  as  shown  by 
the  solid  line  in  the  following  Figures  25  through  36. 
Comparisons  between  experiment  and  numerical  (Appendix  E) 
results  are  presented  for  total  pressure  ratio,  Mach  number 
and  flow  angle.  Figures  25  through  28  show  the  total  pressure 
ratio  comparison.  It  is  noted  that  uncertainty  in  the  probe 
measurement  of  total  pressure  is  small.  The  probe  calibration 
uncertainty  enters  only  when  comparing  Mach  number.  Wake- 
profile  phase  coincidence  was  established  using  the  50  percent 
span  Mach  number  data.  The  phase  was  kept  constant  for  the 
circumferential  comparisons  at  all  other  spanwise  locations. 
Mach  nximber  (Appendix  B)  comparisons  are  contained  in  Figures 
29  through  ^2,  and  Figures  33  through  36  contain  the  flow 
angle  data. 

CFD  total  pressure  ratio  predictions  were  within  10 
percent  of  the  experimental  values  at  25  percent  span.  The 
difference  decreased  to  below  one  percent  at  90  percent  span. 
Mach  number  comparisons  with  measurements  were  within  15 
percent  at  25  percent  span  and  decreased  to  three  percent  at 
90  percent  span.  Flow  angle  comparisons  were  within  five 
degrees  at  90  percent  span  with  improvement  to  within  two 
degrees  at  23  percent  span.  The  reason  for  closer  agreement 


in  prMSur#  rat^io  and  Mach  nuabar  toward  tha  tip  and 

ii^roving  agraaaant  In  flow  angla  toward  the  hub  is 
unaxplainad. 
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rigur«  25.  Total  Pressure  Ratio  Comparison,  25% 
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Figure  26.  Total  Pressure  Ratio  Comparison,  50% 
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Total  Pressure  Ratio 

Corr^sarison  (75%  Span) 


Vigura  28.  Total  Pressure  Ratio  Comparison,  90% 
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MachNuntber 


Wake  Flow  Angle 
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rigurn  33*  Flow  Angle  Conparison,  25% 
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Figure  36.  Flow  Angle  Comparison,  90% 


▼.  C0MCLU8Z0M8  MID  RBC01I1IBMDATZ0N8 

Bxp«riBental  neasurenents  and  computational  output 
locations  were  made  to  coincide.  Repositioning  the 
measurement  voltime  to  one-half  axial  chord  downstream  from  the 
blade  trailing  edges  improved  the  computational  grid 
resolution.  Flow  solutions  with  the  revised  grid  more 
realistically  resolved  endwall  boundary  layer  characteristics 
so  that  future  endwall  LDV  measurements  will  have  a 
corresponding  computed  solution  with  equivalent  resolution. 

Pressure  probe  surveys  will  continue  to  be  important  in 
comparing  experiment  with  CFD  prediction.  The  cobra  probe 
calibration  obtained  was  reliable  and  repeatable.  An 
extensive  set  of  baseline  data  has  been  accumulated  for  future 
use.  The  Scanivalve  data  acquisition  system  demonstrated  the 
ability  to  obtain  pressure  readings  with  a  standard  deviation 
of  less  than  0.2  percent.  The  HP-9000  output  was  manually 
entered  into  a  spreadsheet  program  for  post  processing.  The 
time  required  to  manually  enter  the  output  was  overshadowed  by 
the  spreadsheet  data  manipulation  flexibility. 

The  most  definitive  comparison  between  CFD  and  probe 
measurements  was  in  the  total  pressure  ratio .  Probe 
calibration  error  was  not  present  in  this  comparison.  The 
total  press  are  ratio  disagreed  by  ten  percent  at  the  25 
percent  spar  location  and  less  than  one  percent  at  the  90 
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p^rcMit  location.  In  all  casaa,  tha  coaputer  solution 
pradictad  a  lowar  ratio.  Tha  diffaranca  aay  ba  attributed  to 
conputational  modal ing  paramatars  or  bias  introduced  in  the 
experimental  pressure  ratio  setting.  Upstream  total  pressure 
was  measured  at  only  one  station  and  this  value  may  not  have 
been  a  proper  reflection  of  tha  avarage  total  pressure.  Mach 
number  comparison  characteristics  followed  the  trend  of  the 
total  pressui'e  ratio  as  expected.  Flow  angle  comparison  with 
CFO  was  best  at  25  percent  span  with  a  maximum  difference  of 
5*  at  the  9C  percent  span  location.  Of  all  the  measurement 
locations,  the  90  percent  span  location  was  the  best  defined 
and  flow  angle  measurements  made  there  were  repeatable  within 
0.2  degrees. 

Endwall  LDV  through  a  small  (1.0922  mm)  casing  access  hole 
has  been  demonstrated.  The  ability  to  resolve  endwall 
boundary  layers  and  eventually  tip  flow  losses  appears  to  be 
promising.  Precise  LDV  probe-volume  positioning  techniques 
were  not  obtained.  The  maximum  pressure  ratio  attained  for 
endwall  LDV  was  0.974.  CFD  and  pressure  probe  data  were 
accximulated  at  an  operating  pressure  ratio  of  0.68.  Low  data 
rate,  vibration  of  the  endwall  hole,  and  thermal  expansion 
problems  will  have  to  be  dealt  with  as  the  operating  pressure 
ratio  is  decreased  to  0.68. 

Future  efforts  with  the  computational  fluid  dynamics 
should  include  modeling  the  inlet  bellmouth  and  center  body. 
The  parameters  selected  for  use  in  RVC3D  should  also  be 
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•xplorttd  to  obtain  a  qualitative  feel  for  cause  and  effect  on 
the  solution  with  particular  attention  to  the  unsteady  region 
of  the  exit  flow. 

Pressure  data  acquisition  can  be  continued  to  obtain  a 
nore  conplete  nap  of  the  wake  region.  Pressure 
instrumentation  of  a  blade  and  passage  would  provide  an 
additional  source  of  CFD  code  validation  data.  Data  reduction 
and  presentation  would  be  Improved  with  a  change  in  operating 
system.  The  current  Scanlvalve  system  with  the  HP-9000  does 
not  allow  file  transfer  to  other  computers  for  post 
processing.  HP-9000  limitations  restrict  output  flexibility. 

The  annular  turbine  cascade  now  has  off-angle  LDV  options 
that  have  not  been  explored.  Three-dimensional  radial  LDV 
surveys  of  the  wake  passage  should  be  attempted  for  comparison 
with  pressure  measurements  already  obtained.  Initial 
indications  show  that  the  LDV  probe-volume  energy  may  not  be 
able  to  provide  return  signals  of  sufficient  intensity  as  the 
flow  velocity  increases.  The  addition  of  corrective  optics  to 
the  fiber-optic  probe  may  be  required.^ 

The  endwall  technique  developed  in  this  study  did  not 
provide  pressure  equalization  across  the  access  hole.  The 
pressure  across  the  access  hole  should  eventually  be  equalized 
to  eliminate  flow  through  the  hole  and  therefore  any  influence 
the  opening  might  have  on  the  flow  structure  [Ref.  9].  The 
hole  acted  as  a  natural  optical  noise  filter  and  further 
investigation  of  hole  diameter  Influence  on  LDV  performance  is 


varrant«d.  In  order  to  truly  understand  the  hols  influence 
and  BeasureBtent  reliability,  the  glycerin  seed  particle 
distribution  and  a  study  of  their  aomentuiB  characteristics 
should  be  undertaken,  with  these  refinements,  access  to  the 
flow  through  a  hole  in  the  casing  is  feasible,  without 
changing  the  endwall  flow  characteristics. 
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MTPBMDZZ  B.  COBRA  BROBB  CALZBRfTZOH 

The  following  synopsis  was  extracted  from  Appendix  A  of 
Reference  8.  Mach  number  or  dimensionless  velocity  can  be 
related  to  a  non*-dlmenslonal  pressure  coefficient,  as 
follows: 


Total  velocity: 


Dimensionless  velocity: 


Dimensionless  velocity  and  Mach  number  are  related  as  follows; 


CpT^-CpT 

CpTt 


x.mEL 

v'2^ 

so  that 

- 2 - 

2 


or 


2X^ 

-1)  (1-^2) 
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H  pnMMttTtt  eo«fflei«nt  for  thm  probo  aligiMKl  in  yaw  can  ba 
written  aa 

a  j£-:£sYfiY££V^\ 

P-  APeJVPAYM*) 

where  and  p„  are  probe  sensor  measurenents .  Using  p,  and  p^s 
respectively,  it  can  be  shown  that 

P^EiIhlmf{M,  Pitch) 

Pi 

in  the  range  of  velocities  which  are  of  interest. 

Figure  Cl  contains  the  calibration  data  without  curve 
fitting  and  Figure  C2  shows  the  validation  data  plotted 
against  the  sixth  order  and  linear  curve  fits  at  zero  pitch 
angle.  The  sixth  order  curve  fit  equation  is 

X=(aO+aip+a2P*+a5P^+a4P^+a5P*+a5P*)xlOOOOO 
with  the  coefficients  shown  in  Figure  Bl. 
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Figure  Bl.  Cobra  Probe  Calibration  Coefficients 
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Oftta  «oquislti<Mi  vas  aecoaplishad  using  an  HP-'9000 
ooaaputar  systaa  with  tha  following  prograa: 


10  !FIIE  NAME  :  COGR^XTC 

20  IDISK  lABEL:  "/AUSTIN^ 

30  1  UT  JEFF  AUSTIN 

40  HAST  NDDEFIED  02/03/94  EV  UXR  JOE  SFTTZ 

50  1 

60  1  IKIS  HOGRAN  RBOQRDS  AND  REDUCES  MEASURED  IRBSSURES 

70  !  IRCN  A  scanivalve  OQMNEL'ig)  TO  A  OHBEE  HOIE  IRESSURB 

80  !  nOK.  Scanivalve  FCRT  OCIWBCTIGNS  ARE  AS  FQLLCNS: 

90  1 

100  !  IRXRMf  \SffiIABC£S  AND  FORT  (XmEdlCNS. . .  (#3  not  used) 

110  !  V  «  DESIRED  S/V 

120  1  A  «  rssbh  s/v  fort 

140  1  C  «  CHANNEL 

150  1  VO  -  EKM  READING 

160  1  S$  »  SCANNER  LTSTEN  CODE 

170  1  S  »  SCANNER  # 

180  !  Scanivalve  #5  SET  TO  A1A;ANCE  GNL^ 

190  1  Banb  B  Anfeient  I^nessure  (psi) 

200  i  Patm  B  Atinoepheric  pressure. ..Port  #  1 

210  i  Peal  B  Calibration  pressure. ..Port  #  2 

220  ]  Ptot  B  iteference  total  pressure. .  .Port  #  4 

230  1  Rnjib  *  mner  hub  static  pressure. .  .Port  #  5 

240  1  Pstat  B  TJpstxeam  static  pressure. .  .Port  #6 

250  1  PI  B  Center  hole  on  probe... Port  #  7 

260  1  F23  B  Yaw  port  pressure. .  .Port  #  8 

270  !  beta  »  (P1-P23)/P1 

271  !  prat  B  Riub/Pt^  (pressure  ratio) 

272  !  mach  »  nach(Pl,Fhub,GaiiiDa) 

280  i 

290  OPnON  BASE  1 
300  GaninaFsl.4 

310  CLEAR  Sa^EEN 

311  IKENISR  IS  CRT 

320  DISP  •*P[£ASE  WAIT  NHIIE  RESETTINS  ScanivalveS" 

330  HONT 
340  V*5 
350  AfI 
360  G06UB  Read 

380  IRINT  "DATA  AOQUISITZON  FOR  THREE  HOIE  YAH  IROEE" 

390  DIM  Y(20) 

400  DIM  Patai(10,10),  Pcal(10,10)rPtot(10,10),mib(10,10),Pstat(10,10) 
410  DIM  Pl(10,10),  F23(10,10) 

420  1 

430  INHJT  "IMERICNIH,  DAY,  YEAR  (IE.  02,05,94)"  ,Y(3) ,  Y(5) ,  Y(7) 
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431  INHJT  "INEER  BUN  #:  •*,  Y(9) 

432  FRnnSl  IS  711 

440  mnrr  using  "k,cd,"V''"»dd» of  hun:'',y(3),y(5),y(7) 

441  RONT  "lAm  FUN  ",Y(9) 

450  PRINT 

451  FRnnSR  IS  CRT 

460  !***************************************************************** 
550  HUNT 

560  PRINT  "CCMSTANT  FLON.  WAKE:  49  TO  CGGREES.  SPAN:  5  lOCATIONS 

570  PRINT 

571  PRINT  ••  Zero  EWM  on  Scanivalve  #5,  port  #1.” 

572  PRINT 

580  PRINT  "  Set  20  inHg  for  ScEOiivalve  czdibration.” 

581  PRINT 

582  PRINT"  Adjust  ^>an  for  CKM  »  0.009823  volts  on  port  #2." 

600  PRINT 

610  PRINT"  INSURE  lAM  ADJUSTED  AND  EESIRED  PRESSURE  RATIO  SET. . ." 

630  FOR  KAKE  «  1  TO  10 

640  INPUT  "INTER  AMBUNT  PRESSURE  IN  ir)Hg",Plnhg 

641  P^aiEb=Pinhg*(.  4911541) 

650  PRINT 

660  PRINT  "CHECK  DESIRED  SPAN  AND  FLON  ALKaOIENT  SET. . ." 

661  PRINT 

662  PRINT  «*****FRESS  OONTINUE  WHEN  READY  TO  TAKE  DATA*****" 

670  PAUSE 

671  CLEAR  SCREEN 

680  PRINT 

700  FOR  SPAN  =  1  TO  10 

710  V=5 

720  J^l 

730  G06U3  Read 

750  ! 

850  IRBOORD  THE  PRESSURE  DATA  FOR  THIS  WAKE  AND  SPAN 
860  ! 

870  ! 

880 
890 
900 
910  ! 

912 

913 

914 

915 

916 

919 

920 

921  ! 

922  ! 


OUTEUT  722;"F1R1T1M0Z1  "  !  SETS  UP  HP3456A  DVM 

!  Scanivalve  #5  DESIGNATED 
S=1  !  HP3495A  SCANNER  #1  DESIGNATED 

PRINTER  IS  711 
PRINT 

PRINT  "Wake  Data  Groip  ",  WAKE 
PRINT  "Span  Location  Ocunter  ",  SPAN 
PRINT  USING  "/,5X,K,DDD,/,2X,K,15X,K,17X,K";" 
Scanivalve  #",V,"It>rt", "Volts", "Psi" 

PRINT 

PRINTER  IS  CRT 
RE^JD  ScanivcQve  PORTS 
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923  1 

930 

940 

950 

960 

970 

980 

990 

1000 

1002 

1010 

1020 

1021 

1022 

1023 

1024 

1025 

1026 
1027 

1030 

1031 

1032 

1033 

1034 

1035 

1036 

1037 

1040 

1041 

1042 

1043 

1044 

1045 

1046 

1047 

1050 

1051 

1052 

1053 

1054 

1055 

1056 


VCR  Ar-1  TO  8 
OOSUB  RmkI 
IflOT  2.0 

anivr  7oi;''c"  i  clears  scanner  «i 

cxmur  701  using  •myrVfS  l  scanner  set  to  read  scanlvalve 

TRIGGER  722 
bhqi  722  ;vo 

HmnER  IS  711 
!  Pressure  in  psi 
IF  AfI  them 

Pfttn(1AXB,SPAN)^V0*1000 

ERUir  USING  "2X,nX>,7X,7DD.I]0CmX>,8X,7D0.DIXD''; 

P,  VO ,  pBt3B(NAKE,  SPAN) 

HUNIfR  IS  CRT 

PRINT  USING  "2X,IXX3,7X,7DD.DC)DDDDD,8X,7DD.DDC)D''; 

P,V0,Pat3n(NAXE,SEAN) 

ERINTER  IS  711 
END  IF 
I 

IF  Af2  then 

(lOKE,  SPAN)  aVOnOOO 

PRIHr  USING  "2X,I]DD,7X,7I».DDDDDDD,8X,7DD.DDDD''; 

P,  VO ,  Peal  (HAKE ,  SPAN) 

FRnnER  IS  CRT 

IRIMr  USING  "2X,DIX5,7X,7DD.im»DD,8X,7DD.DDDD’’; 

P,VO,Pcal  (NAKE,SPAN) 

HRUflER  IS  711 
END  IF 

! 

IF  Ai=4  THEN 

Pbot  (HAKE,  SEAN)  -(VD*1000)  -Patin  (WAKE,  SPAN)  +Pairb 
ERINT  USING  "2X,IXX>,7X,7n>.im]DClD,8X,7DD.DIX)D’'; 

P,VO,Ptot(WAKE,SPAN) 

IRINIER  IS  CRT 

ERIN!  USING  ''2X,C0D,7X,7DD.I]E)DDDDD,8X,7DD.DDDD"; 

P,VO,Ptot(WAKE,SPAN) 

IRINZER  IS  711 
END  IF 
I 

IF  Af5  THEN 

Anib  (WAKE,  SPAN)  »(V0«1000)  -Patm  (WAKE,  SPAN)  +Panto 
IRINr  USING  ''2X,IXO,7X,7DD.DDDDDDD,8X,7DD.DDDD''; 

P,  VO ,  ItlUb  (WAKE,  SPAN) 

IRINIER  IS  CRT 

ERINT  USING  "2X,DC»,7X,7DD.DDim)D,8X,7DD.DDC)D''; 

P,VO,Itlub(WAKE,SPAN) 

IRINIER  IS  711 
END  IF 
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1057 

! 

1060 

IF  TOEN 

1061 

PBtat  (fAKE,  SFAN)«(V0*1000)  -Fatm  (WAKE,  SPAN)  +Panb 

1062 

Rmfr  USING  "2X,CaD,7X,7IX>.D0aDCa3,8X,7IX).DDt]D”; 

P,  VO,  PBtat  (HAKE,  SPAN) 

1063 

FKQfllK  IS  CRT 

1064 

HOOT  USING  '*2X,DDD,7X,7DD.DDDDDtX),8X,7DD.DDDD"; 

P,  VO ,  PBtat  (HAKE,  SPAN) 

1065 

HUNIfR  IS  711 

1066 

QID  IF 

1067 

1 

1070 

IF  A=7  IHEN 

1071 

PI  (HAKE,  SPAN)  «(V0*1000)  -PBtm(HAKE,  SPAN)  +Pairi3 

1072 

PRINT  USING  "2X,CGD,7X,7EO.IXaXXD,8X,70D.DC0D"; 
P,V0,P1(HAKE,SPAN) 

1073 

IRUaXR  IS  CRT 

1074 

PRINT  USING  ••2X,Dl»,7X,7I».DDra3DD,8X,7DD.DDDD''; 
P,V0,P1(WAKE,SRAN) 

1075 

PRINTER  IS  711 

1076 

END  IF 

1077 

1 

1080 

IF  Af8  THER 

1081 

P23  (HAKE,  SPAN) = (V0*1000)  >Fatin  (HAKE,  SPAN)  +Pairb 

1082 

PRINT  USING  ''2X,DDD,7X,7DD.DCCaDCD,8X,7DD.DC0D"; 

P,  VO,  P23  (HAKE,SPAN) 

1083 

PRINIXR  IS  CRT 

1084 

PRIWr  USING  •'2X,Dl»,7X,7DD.DDDDDDD,8X,7DD.DDDD''; 

P,  VO ,  P23  (HAKE,  SPAN) 

1085 

PRINTER  IS  711 

1086 

1 

1087 

prat?>Ptiub(HAKE,SPAN)  /Ptot  (HAKE,  SPAN) 

1088 

beta=  (PI  (WAKE,  SPAN) -P23  (WAKE,  SPAN) ) /PI  (WAKE,  SPAN) 

1089 

inach»SQRr(5* (Pl(HAKE,SPAN)/Phub(HAKE,SPAN) )  ^  (1/3 . 5) 

1090 

END  IF 

1091 

1 

• 

1092 

NEXT  A 

1097 

PRENT  "prat  =  ",  prat 

1098 

PRINT  "beta  *=  ",  beta 

1099 

PRINT  "mach  §  -  ”,  roach 

1100 

PRINIER  IS  CRT 

1101 

PRINT 

1102 

PRINT  "Set  Span  and  align  prcbe. .  .CCXniNUE  WHEN  READY. 

1103 

PRINT 

1104 

PAUSE 

1105 

CLEAR  SCREEN 

1106 

NEXT  SPAN 

1107 

PRINT  "Verify  wate  &  vpdatB  Panto. .  .CONTINUE  WHEN  READY. 

1108 

PRINT 

1109 

PAUSE 

1110 

CLEAR  SCREEN 

1111 

NEXT  HAKE 
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********«********B4D  OF  MAIN  PBXaWt*********************** 


m2 

1113 

m4 

U15  I  SUERXJnNE  TO  FQSmCXI  AND  READ  S/V  PORTS 
1116  1  Scanlvalvft  il  HAS  IHOP  HOME  FEATURE 

m?  RaadraTnvr  707  using  '<,k";v 

1118  P0i«P0U.(707) 

ms  L4CIMAND(F0,15) 

1120  Tl-imXIE(P0,4) 

im  M1-BIMAMD(T1,7) 

1122  M0«M1+L 

1123  CLEAR  707 

m4  IF  THEN  Finish 

1125  GUnUT  701;  “C" 

1126  OmUT  701  USING  "DD»‘;V-1 

1127  OUmiT  701;  "C" 

1128  MATT  .1 

1129  GOTO  Reed 

1130  FinishzREIURN 

1131  END 


APPBMDZZ  C.  ORZD  QEXZlUkTZON 


Grid  generation  was  conpleted  on  IRIS  workstations. 
TCGRID,  a  FORTRAN  progran,  required  the  adjustment  of 
parameter  statements  to  dimension  arrays  for  the  model  build. 
Parameter  adjustment  required  patience.  The  author  found  the 
manual  helpful  in  most  cases,  however,  the  most  useful  tool 
turned  out  to  be  a  running  summary  of  the  namelist  input  for 
troubleshooting.  The  most  difficulty  was  observed  in 
obtaining  a  smooth  solution  on  the  suction  side  of  the  blade 
geometry.  The  following  TCGRID  namelist  file  was  utilized  to 
obtain  the  grid  shown  in  Figures  Cl  and  C2.  Reference  4 
contains  the  original  list  of  blade  coordinates  and  is  not 
included  here. 

&naml  iiip=150  jm=31  km=65  itl-20  ics^lS  k2(^3  merid=0  &end 
&nam2  nl^l6  nteF^io  dsl^.018  dsteF=.003  dshubB.0004  dstif^.0004 
dsiitep^.OOl  dswex^.060  dsthz^l.  dsndm.0004  dsraax=.025 
dsra=.45  rcon^.OSS  6&nA 

&nam3  itetn»100  idbgpO  0  0  0  0  0  0  0  aabb=.5  codd=.45  &end 
&nam4  zbo»-1.5  -1.5  0.5  -1.5  -1.5  .5 

zbo«  3.895  3.895  3.895  4.585  4.585  4.585  &end 
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Pigura  C2.  Grid  Geometry,  3-D  Perspective 
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APFBMDZZ  D.  KVCSD  OODI  ZMFOT 


RVC3D  runs  were  initially  coaplstsd  on  IRIS  workstations 

at  NFS.  Solutions  were  United  to  less  than  500  iterations  to 

ensure  residuals  were  decreasing.  File  transfer  of  the  fort.l 

(grid)  and  the  rvc3d.in  (naaelist)  files  to  the  Cray  were 

completed  only  after  Initial  testing  was  completed.  The 

following  namelist  contains  the  commands  used  for  the  final 

set  of  iterations  to  the  minimum  residual  attainable.  It  is 

a  restart  (iresti=l)  from  1200  iterations  with  termination  at 

1490  iterations  (imax=290). 

'IRANSGNIC  lURBINE  Annular  Cascade  >  (Restart  12>1490) ' 

&nll  ijv«150  jiB=31  kiff=€5  itl»20  iil»67  fiend 

finl2  cl»=5.5  aviscl=0.0  avisc2=0.0  aviso4=0.50  ivdt?*l 

nstig=4  itmax>290  Itb^I  epi=°O.50  epjsO.60  eplc=0.60  fiend 
finl3  IbclnFl  iboexB3  isymb=0  ireesio  icmteSO 
iresti»l  iresto^l  igin=0  fiend 

finl4  eni»e=0.13  enity=0.0  enirzM}.0  e9qpb=0.0  prab=0.6800 
ga^l.4  an=0. 000000  igecnt=l  2de>p»-67.0  fiend 
finl5  ilt?=3  tip^l.OO  rair=6.651e6  pmr^.7  prtr=.9 

vispwTB. 666666  srtip=0.0  aiut]D»14.  jedge^31  kedge=20 
iltiir*2  dblhFO.0024  dblt=0.0089  fiend 
finl6  iol=l  io2»165  oar>0.  ixjl^  njo=l  hkD=3  jo=l 
ko»5  11  16  fiend 

Runs  on  the  NFS  cray  were  accomplished  with  the  command 
>gsub  myjob.ngs 

where  the  file  myjob.ngs  was  as  follows: 
myjob.ngs 

#QSUB  -g  reg  -IT  43200  -IM  9MW 
cd  cfd 

rvc3d  <rvc3d.in>  rvc3d.out 
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Th«  RVC3D  axecutable,  grid  (fort.l),  restart  (fort. 2), 
and  the  Input  nanelist  (rvc3d.in)  were  held  in  the 
subdirectory  cfd.  The  final  solution  required  13.2  hours  of 
Cray  computer  time.  A  plot  of  the  residuals  is  shown  in 
Figure  Dl. 


RVC3D  Residuals 

Normalized  to  Maximum 


Figure  Dl.  RVC3D  Residuals 
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JkVPBMDXX  B.  B¥C3D  DBTA  BBDUCTZOM 


Exit  plane  data  were  reduced  using  the  FORTRAN  program 
plane. f.  RVC3D  output  the  following  non-dimensional  values 
for  each  grid  point  in  the  file  fort. 3: 


Qi 


=-£. 

Po 


Q2  = 


9oCq 


<73 


Qa 


Po^o 

_  pw 


Po^o 


05  =  -^ 

Po 


Co’  c. 


Computational  output  was  manipulated  in  plane. f  for 
comparison  with  experimental  results.  The  following  table 
provides  a  coordinate  system  comparison. 


Plot3d 

Cartesian 

Annular  Cascade 

Cartesian 

cylindrical 

x(i)  -  u 

z(i)  -  q2 

z(i)  -  q2 

y(i)  -  V 

y(i)  -  q3 

V«(i) 

x{i)  -  04 

V.(i) 
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Transformation  to  cylindrical  coordinates  was  conducted 
using  the  grid  location  information  contained  in  the  RVC3D 
grid  input  file  fort.l. 

6{i)  -sin~^  ]  r{i)WynP^^TIP 

All  angles  (6)  were  referenced  to  the  y  axis  shown  in  the 
blade  geometry  plot  of  Figure  3.  Flow  angle  computations  were 
completed  using  the  non-dimensional  q  values  described  above 
as  follows: 

Flow  Angle  T  (dimensional) : 

Ve=- (u)  sin0+ (v)  COS0  V^=  (u)  cos0+ (v)  sin0  V,=w 

Flow  Angle 


Direct  comparison  was  possible  with  pressure  ratio 
measurements.  The  ratio  between  local,  and  upstream  reference 
total  pressures  was  computed  as  follows: 


Po 


=Y(Y-1) 


2  Po 
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■  JE,  ♦  .»i_  p(u  * + V*  ♦  V*' 

Po  Po 
^•Y  (Y-1) 


Mach  ntiaber  was  obtained  using  the  isentropic 
relationship  between  Mach  number  and  temperature  ratio. 


g  =  — — _  _?! + — +g? +g4 ) 

y(y-i)  Tq 


The  program  plane. f  is  included  here  for  reference.  It 
is  an  updated  version  of  the  program  written  by  Thomas. [Ref. 
4  ] .  The  new  version  includes  a  calculation  of  flow  angles 
after  conversion  to  cylindrical  coordinates  for  proper 
comparison  with  experimental  results. 


c********************************************************************* 

c 

c  plane,  f  reads  xvc3d  files  &  writes  ascii  files  for  plotting 

c  exit  plane  mach  nusber  and  flow  angles 
c  (ccgaparison  with  eigjerimented  results) 

c  unit  1  «  input  xyz  file 

c  unit  3  »  input  q  file 

c  unit  4  s  output  residuad  history 

c  unit  7  «  output  (k  #1)  counter,  theta,  pt/pO,  Mach,  Angle 


c  unit  8  «  cutput  (X  #2)  ocuntar,  theta,  iX/pO#  Mach,  Angle 

c  unit  9  ■■  output  (k  #3)  oounter,  theta,  pi/pQ,  Mach,  Angle 

c  vnit  10  -  output  (k  #4)  oounter,  theta,  pt^^,  Mach,  Angle 

c  uilt  11  «  a^tput  (k  #5)  oounter,  theta,  ptv^,  Mach,  Angle 

o 

pazameter  (ni«150 ,  njB3 1 ,  rik»65 ) 
integer  vm 

reed  x(ni,ir),hk),y(nl,nj,rik),z(ni,nj,rik),  gama 
reed  qg(5,nl,nj,hk,read(5000,5)  ,gi,q[2,(^,q4,g5,in(5) 
dimension  ](k(5)  ,ang(5)  ,theta(5)  ,tt0(5)  ,vtnd(5)  ,vznd(5} , 

1  V2na(5),ptp0(5),deg(5) 
o 

o  k-^values  are  harcHidred  below  (10%,  25%,  50%,  75%,  90%) 
o  (use  i^)an.f  to  obtain  %-k  equivalency) 

o 

data  k.V21,27,33,39,45/ 

c********************************************************************* 

o  read  grid  ooordinates 

c********************************************************************* 

read  (1 ,  *)  im,  jm,  km 

read(l,*)  (((x(i,  j,k)  ,i*l,im)  ,j*l,jm)  ,Jc=l,lan) , 

1  (((y(i,j,k)  ,i=l,im)  ,j=l,jm)  ,k«l,km) , 

2  (((z(i,  j,k)  ,i»l,im) ,  j=l,jm)  ,k=l,km) 

c********************************************************************* 

o  reed  restart  file 

c********************************************************************* 

read  (3 ,  *)  imax,  jinax,kmax 
read(3,  *)  fsnach,alpha,re,tiiae 
o 

icheck^iabs  ( im-imax)  +iabs  ( jm-jinax)  +iab6  (km-kmax) 
if  (icheck.ne.  0)  then 
write(6,610)  im,  jm,km,imax,jmax,kinax 
stop 
endif 
c 

read(3,*)  ((((qq(l,i,  j,k)  ,i=l,im)  ,j=l,jm)  ,k=l, km)  ,1=1,5) 
o 

o  additioned  residued  data 

read(3,  *)  itl,  iil,phdeg,gama,cni,nres,dum,dum,dum,dum 
reai^(3»*)  ((i«sd(nr,l)  ,nz=l,nres)  ,1=1,5) 

ctt******************************************************************** 

c  Pressure  side  of  exit  wake  cut. . . 

l^kk(3) 

i=l 

do  10  j=jmax,l,-l 
UDF^mH-l 
do  15  1=1,5 
kf=kk(l) 


c 


70 


C 

C 


c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 


nariiwa  daslred  q  values  for  use 

qi-qq|(l,i,j,k) 

q2-cn(24»j«K) 

q3-qq[(3,l,j.]c) 

q4-qq(4,i»j/)c) 

q5i^(54#j,K) 

Obspute  pt/pO  ratio 

ptpO  (1)  xqana*  (geuna-l)  *q5+(gaflBa-gaaB**2/2)/qi*  (q2**24q3**24q4**2) 

Ocepite  T/Tt  ratio  and  nach  twber 

ttO  ( 1 )  -^paa*  (gama-l)  /qi*  (q5-0 . 5*  (q2**2-»q3**24q4  **2 )  /ql) 

»(!)-( (l/ttO(l)-l)  *2.0/ (gana-l)  )**0.5 

Ccnfxite  fOieta  zmd  Non-D  velocities  theta,  r  and  z 

theta(l)=asln(y(i,j,k)/(z{i,j,k)**2+y(i,j,k)**2)**0.5) 
vtnd(l)«-q4*sin(theta(l)  )4^*oo8(theta(l) ) 
vmd  (1) ’>^*006 (theta  (1)  )4q3*sin  (theta  (1) ) 
vznd(l)<>q2 

deg(l)*theta(l)  *57.29f 

Oonpute  flow  angle  (degrees)  using  vtnd  and  vznd 

ang(l)»atan2(vtncl(l)  ,vzna(l))*57.296 
15  oontinue 

write(7,300)um,deg(l)  ,ptpO(l)  ,n(l)  ,ang(l) 
write(8,300)um,deg(2)  ,ptpO(2)  ,in(2)  ,ang(2) 
write(9,300)uni,deg(3)  ,ptp0(3)  ,n(3)  ,ang(3) 
write(10,300)um,deg(4)  ,ptpO(4}  ,in(4)  ,ang(4) 
write  (11, 300)  urn,  deg  (5)  ,ptp0(5}  ,]n(5)  ,ang(5) 

10  oocitinie 

Repeat  above  on  other  side  of  wake  cut  at  inax 
jsl  not  used  (di^licate  location  across  wake) 
jnax  not  repeated  (same  location  as  j=31  at  i=l) 

jinaxnr=^]nax-l 

i«=ijn 

do  30  j»2,  jniaxm,! 
unF^nH-l 
do  25  1=1,5 
k=kk(l) 

ql=qq(l,i,j,k) 

^=qg(2,i,j,k) 

‘i3=*}i(3,i,j,k) 

q4=qq(4,i,j,k) 

q5=q[j(5,i,j,k) 
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ptpO(l)>^|aina*  (gana-l)  *<J5+(ganB-gaitB**2/2)/ql*(<p**24<j3**2-K34**2) 
ttO  (1)  •‘Cfama*  (geana-l)/ql*  (q5-0. 5*  (cp**2+^**24q4**2)/ql) 

»(!)-( ( VttO(l)  -1)  *2/  (gasa-l) )  **0.5 
theta(l)-asln(y(l,j,k)/(z{i,j,k)**24y(i,j,k)**2)**0.5) 
vtnd(l)»-q4*sin(theta(l)  )4q3*aoB(th0ta(l) ) 
vmi(l)»q4*oos(the!ta(l)  )4q3*sin(theta(l) ) 
vznd(l)^|2 

dea  ( 1)  "theta  ( 1)  *57 . 296 
ang(l)"atan2(vtnd(l)  ,vznd(l))  *57.296 
25  cxxitinLie 

write(7,300)ian,deg(l)  ,ptpO(l)  ,in(l)  ,ang(l) 
wite(8,300)uni,deg(2)  ,ptp0(2)  ,n(2)  ,ang(2) 
write(9,300)um,deg(3)  ,ptp0(3)  ,ai(3)  ,ang(3) 
write  (10, 300)  urn, deg(4)  ,ptp0(4)  ,in(4)  ,ang(4) 
write  (11, 300)  deg  (5)  ,p4)0(5)  ,b(5)  ,ang(5) 

30  oontinue 

c  residual  history  out^xit  to  unit  4 

C*A**Ait************AAA*it**4t**4t  **************************************** 

write(4,310)  1, (resd(l,l),l=l,5) 

do  40  j»2,nres 

it=l0*(j-l) 

40  write(4,310)it, (resd(j,l),l=l,5) 

c*************************.^******************************************* 
300  format(i5,5f8.3) 

310  fonnat(i5,5(lx,el0.3)) 

610  fonnat('  *****  warning  ******,/, 

1  '  im,  jm,  km,  read  frcro  input', 3i5, ^  do  not  match',/, 

2  '  im,  jm,  km,  read  frcni  restart  file',3i5) 
stop 

end 
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